import numpy as np
import math

L=4
jiange=20
tota=np.loadtxt("E.dat")
for ig in range(100):
    ihmin=ig*jiange
    ihmax=(ig+1)*jiange-1
    a=tota[ihmin:ihmax,1]
    b=np.mean(a)
    sd=0
    n=len(a)
    for i in a:
        c=(i-b)**2
        sd=sd+c
    f=open('sigmaE'+str(L)+'.dat',mode='a+')
    f.write(str(tota[ihmin][0])+'\t')
    f.write(str(b)+'\t')
    f.write(str(np.sqrt(sd)/np.sqrt(n*(n-1)))+'\n')
    f.close()
    
    
tota=np.loadtxt("CV.dat")
for ig in range(100):
    ihmin=ig*jiange
    ihmax=(ig+1)*jiange-1
    a=tota[ihmin:ihmax,1]
    b=np.mean(a)
    sd=0
    n=len(a)
    for i in a:
        c=(i-b)**2
        sd=sd+c
    f=open('sigmaCV'+str(L)+'.dat',mode='a+')
    f.write(str(tota[ihmin][0])+'\t')
    f.write(str(b)+'\t')
    f.write(str(np.sqrt(sd)/np.sqrt(n*(n-1)))+'\n')
    f.close()


tota=np.loadtxt("M.dat")
for ig in range(100):
    ihmin=ig*jiange
    ihmax=(ig+1)*jiange-1
    a=tota[ihmin:ihmax,1]
    b=np.mean(a)
    sd=0
    n=len(a)
    for i in a:
        c=(i-b)**2
        sd=sd+c
    f=open('sigmaM'+str(L)+'.dat',mode='a+')
    f.write(str(tota[ihmin][0])+'\t')
    f.write(str(b)+'\t')
    f.write(str(np.sqrt(sd)/np.sqrt(n*(n-1)))+'\n')
    f.close()
    
    
tota=np.loadtxt("Br.dat")
for ig in range(100):
    ihmin=ig*jiange
    ihmax=(ig+1)*jiange-1
    a=tota[ihmin:ihmax,1]
    b=np.mean(a)
    sd=0
    n=len(a)
    for i in a:
        c=(i-b)**2
        sd=sd+c
    f=open('sigmaBr'+str(L)+'.dat',mode='a+')
    f.write(str(tota[ihmin][0])+'\t')
    f.write(str(b)+'\t')
    f.write(str(np.sqrt(sd)/np.sqrt(n*(n-1)))+'\n')
    f.close()
